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Is the CMB asymmetry due to the kinematic dipole? 
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ABSTRACT 

Parity violation found in the Cosmic Microwave Background (CMB) radia- 
tion is a crucial clue for the non-standard cosmological model or the possible 
contamination of various foreground residuals and/or calibration of the CMB 
data sets. In this paper, we study the directional properties of the CMB parity 
asymmetry by excluding the m = modes in the definition of parity parameters. 
We find that the preferred directions of the parity parameters coincide with the 
CMB kinematic dipole, which implies that the CMB parity asymmetry may be 
connected with the possible contamination of the residual dipole component. We 
also find that such tendency is not only localized at I = 2, 3, but in the extended 
multipole ranges up to I ~ 22. 



Subject headings: cosmic microwave background radiation 
methods: data analysis — methods: statistical 



Introduction 



early universe 



Symmetry of the physical process in our Universe and particular mechanisms of its 
violation is golden mind of the modern physics. Since pioneering Lee and Yang investiga- 
tions of the parity symmetry in the weak interaction, the principle of symmetry is deeply 
incorporated into the modern particle physics, including the Higgs mechanism of symme- 
try breaking, in chemistry, in physics of condensed matter and, in general, in the theory of 
the phase transition. Passing from the microscopic physics to the properties of the space 
and time at large, we have to admit that the Cosmic Microwave Background (CMB) radi- 
ation anisotropy provides invaluable test for the investigation of parity at the megascopic 
scales above the scale of inhomogeneity 



Gruppuso et al.ll201ll ; lHansen et al.ll2011 



(Land k, Maeueiio 


2007; 


Kim & Naselskv 


2010alJbl: 


Maris et al. 


2011; 


Ben-David et al. 


2011). showing 



significant dominance of the power spectrum stored in the odd multipoles over the even ones. 
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Recently, in (IKim &: Naselskyll201ll ) it was shown that the odd multipole preferences tidily 
connected with the anomalies o f the two-point corre l ation function, in particular, the lack of 



correlations at 60° < O < 180° (jSchwarz et al. 



anomaly of the correlations at 1° < < 30° (IKim &: Naselsky 



2004 



Copi et al. 



2010), and newly discovered 
201ll ). In combination with 



the w idely discussed anomalies of the CM B map and the power spectrum ( lEriksen et al. 



20041 ) (see for review (IBennett et al.l l201ll )). investigation of the origin of these anomalies 
could put a new light on the physics of the early Universe, the methods of the foregrounds 
reductions and calibration of the CMB data sets, improving our knowledge of the most 
fundamental cosmological parameters and the theory of inflation. 

The local motion of an observer through the CMB frame produces the so-called Kine- 
matic Dipole (KD) anisotropy of the CMB, which is the most powerful component of the 
signal and fitted out from the CMB data before cosmological analysis. In this paper, we are 
going to show that some of the discovered features of the Wilkinson Microwave Anisotropy 
Probe (WMAP) CMB TT anisotropy, including the odd-parity preference of the power 
spectrum, could have common origin associated with the KD of the CMB. Previously, the 
possible contamination of the CMB by KD has been as sessed by multipole vector statistic 
(align ment of the quadrupole and octupole components) (IGordon et al.ll2005l ; iPeiris fc Smith 
2010l ). Additionally, we will show that the low multipole anomalies are associated with other 
anomalies such as the lack of angular correlation, the even/odd-parity asymmetry and the 
planarity of the multipole I = 5. First, we will focus on the properties of the CMB TT 
correlation function C(0) = AT(n)AT(n') at angle = arccos(n • n') = tt and show that 
C(0 = 7r) is connected to the parity parameter g(l). We will estimate the power spectrum 
without the m = mode so that the rotational invariance of the angular power spectrum 
may not be automatically satisfied. As well known, these m = modes pick up a certain 
direction (i.e. the z- axis direction in the spherical coordinate system where ai m are defined 



(IGordon et al.ll2005l )). Using these estimators, we will investigate the possible preferred di- 
rection in the CMB field. If our Universe is, indeed, statistically homogeneous and isotropic 
with Gaussian seed perturbations, we should have no or little parity preference, given for 
our estimators associated with the angular power spectrum. Using the WMAP 7-yr Inter- 
nal Linear Combination (ILC7) map, we compute the parity parameters for coordinates of 
various orientations. We find some level of alignment between the KD direction and the ori- 
entation, in which the parity asymmetry is greatest. Therefore, the CMB parity asymmetry 
may be related to the systematics associated wit h KD, which may be also responsible for the 



align ment problem of quadrupole and octupole (IGordon et al.ll2005l ; Ide Oliveira-Costa et al. 
2006h . 



The outline of the paper is the following. In Section 2 we introduce the basic char- 
acteristics of the CMB parity asymmetry. In Section 3, we investigate the orientations of 
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maximum parity asymmetry and compare them with the CMB kinematic dipole. In Section 
4, we summarize our investigation. 



2. Odd-multipole preference of the CMB power spectrum 

The temperature fluctuations of CMB anisotropy, can be conveniently decomposed as 
follows: 



oo I 



AT (9, 



yj O'lmXlmi 



1=0 m =-l 



where ai m are the coefficients of decomposition: a/ m = \ai m \ exp(i0; m ), with 0; m as the phase. 
Under the assumption of total Gaussian randomness, as predicted by the large class of 
inflationary models, the amplitudes \ai m \ are distributed according to Rayleigh's probability 
distribution f unction and the ph ases of a/ m are supposed to be evenly distributed in the 
range [0, 2n] ( iBardeen et al.lll986l ). 



For any signals T(n) defined on the sphere, one can extract symmetric (AT + (n) 
AT + (— n)) and antisymmetric (AT~(n) = -AT^(-n)) components, where 



AT 



n 



«^r ± (/)^ m (n) 



(2) 



l m=—l 



and = cos 2 (f ), r-(/) = sin 



2 i nl 



f ), Y lm (h) 



■1YY, 



lm \ 



-n). Naive expectation, where 



the concordant ACDM cosmological model with initial statistically isotropic and Gaussian 
adiabatic perturbations is assumed, is the absence of any features distinct between even and 
odd multipoles. However, in reality this statement needs more accurate clarification. In 
particular, for the curvature perturbations beyond the present horizon the power spectrum 
is given by P( k) oc k~ A+rbB , where n s ~ 0.96 is the spectral index of the primordial density 



perturbations (IKomatsu et al.ll201ll ). Thus, the variance of the metric perturbations a 2 ~ 
J k 2 P(k)dk oc fc^;" 1 has very weak power-law (n s < 1) or logarithmic divergence (n s ~ 1), if 
^min — > 0. Since the low multipole range of the CMB temperature anisotropy is determined 
by the ordinary and integrated Sachs- Wolfe effects, these peculiarity of the power spectrum 
of metric perturbations are crucial for the two-point correlation function: 



C th (8) = < AT(n)AT(n') > 
00 9/4-1 



l=L 



47T 



(3) 
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where Cth{l) are the theoretical power spectrum, P;(cos0) are the Legendre polynomials, 
cosO = n ■ n', and < .. > denotes the average over the statistical ensemble of realizations. 
Using Eq. ([3]), we may easily show, for the largest angular distance: 

00 2/ + 1 

c th (e = n)=Y, -^-cmq+V) - r-(i)). (4) 

' ^min 

As clear in Eq. (151) , the natural way to estimate the relative contribution of even and odd 
multipoles to the correlation function is to define the statistic 

T!v=v . ^c(Z')r+(Z') 
" U Ytv=i>. ^C{V)Y-{l>y (5) 

min ^" 

where l' min = 1 or 2 (see discussion in the forthcoming section s) . Not e that this statistic 
differs from the g(l), widely used in fcim fc Naselskvi boiOaU bl. 120111 ). in the sense that 



l'{l' + 1) is replaced by 21' + 1. Then, from Eqs. (j4]) and (JSJ) we get the estimator of the 
quantity C t h{@ = 7r): 

C(0 = tt) = P-(Z) b(Z) - 1] , 

pHi) = E ^ccoi^CO- (6) 



m 1 n 



Thus, if p(Z) = 1, the corresponding correlation function is C(0 = 7r) = 0. In reality, 
the theoretical correlation function shows some parity asymmetry. For instance, the dipole 
component naturally contributes the odd parity, while the quadrupole contributes the even 
parity. When the background cosmological parameters correspond to the concordant ACDM 
model, the properties of theoretical correlation function C^(9 = tt) depends on the value of 
l' min , which is clearly shown in Fig. [TJ We find that the odd parity is perferred when the odd 
l' min is chosen, while the even parity is perferred when the even l' min is chosen. However, the 
observed data show the different tendency for the parity asymmetry, i.e. the parity violation 
comparing with the theoretical predictions. From Fig. dj we find that the odd parity is 
always perferred for all the cases with l' min < 15. 

Usually, the dipole component of Cth(® = tt) is not included in Eq. and the lower 
limit is set to l' min = 2. In Fig. [2] we carefully calculate the WMAP7 observational correla- 
tion functions by considering the KQ75 mask and a theoretical prediction with la interval, 
where the cosmic variance effect is taken into account. Consistent with Fig. (TJ we expect 
the theoretical value of C t h(0 = tt) is positive, and the zero value is included at 68% C.L. 
Therefore the Universe has tendency to be parity asymmetric g(l) > 1 by choosing l' min = 2, 
with very low chance of g(l) < 1. However, the WMAP7 data show that C{9 = it) < at 



-5 - 



95% C.L., and our Universe belongs to the very rare realizations of the ACDM cosmological 
models with given by the WMAP7 cosmological parameters. 

If the parity violation has the cosmological origin, we should also see the similar parity 
violation in TE, EE and BB components of the polarized signal. However, in the curret 
WMAP data, the noises of the polarization data, including the TE cross-correlation data 
are quite large. So we cannot get any solid results on the parity asymmetry by using the 
WMAP polarization data, when taking into account the large error bars. It is expected 
that the forthcoming Planck data would provide the much better chance to study the parity 
asymmetry in the polarization data, and be helpful to reveal of the origin of the parity 
violation. 



3. Directional statistic of the parity asymmetry 

As shown in the previous section, for the random Gaussian statistically isotropic and 
homogeneous perturbations of the CMB, the correlation function C(O) is fully determined 
by the power spectrum C(l), which is rotationally invariant. Statistical invariance means 
that for any rotations of the reference system of coordinate, the power spectrum and the 
correlation function are invariant. The idea of the method, proposed in this section, is to 
replace the power spectrum C(l) in Eq. (J3J) by a rotationally variant power spectrum D(l), 
defined as 

1 1 

m = 27TT ^ W 2 (i-M, ( ? ) 



where S mm i is the Kroneker symbol. 



As it is seen from the definition given by Eq. ([7j), the relative difference between D(l) 
and C(l) is given by A(Z) = = ~ a m/ J2 m \ a im\ 2 - So we have A(/) ~ O(^) for 

random Gaussian CMB field. Thus, the major difference A (7) comes from / = 2 and I = 3 
modes, while for I > 5 their contributions are smaller than 10%. 

Now, we can study the power spectrum D(l) in any coordinate system. Imagining 
the Galactic coordinate system is rotated by the Euler angle (i/j,9,<f)), and the coefficients 
o; m ('0) 0) 4>) m this new coordinate system can be calculated by 

i 

ai m (ij),0,4>)= "w^L'^^li ( 8 ) 

m'=— I 

where a/ m = a; m (0, 0,0) are the coefficients d efined in the Ga lactic coordinate system, and 



D l mm i(i}),6,(f>) is the Wigner rotation matrix f lEdmondsl Il985l ) . Similar to Eq. (J7|), we can 
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define the power spectrum D(l; if), 8, 0). It is easy to find that D(l; ip, 8, <j)) is independent of 
the angle ip, so in this paper we only consider two Euler angle q = (8, <j)) and set ip — 0. If 
we consider q as a vector, which labels the z-axis direction in the rotated coordinate system, 
then <fi) is the polar coordinate of this direction in the Galactic system 0. 

Now, we can define the rotationally variable parity parameter G(l; q) by replacing C(l) 
in Eq. dSJ) with D(l; q), and estimate the maxima and minima of G(l; q) for different Euler 
angles q. By the definition, the parity parameter G(l; q) depends on the coefficients ajo(q) 
as follows: 

C(l- a) - P+{1) ~ X+{lA) (9) 
where X ± (l;q) = ± Y!v=2 a uo(^) r± ( 1 '), and 

mm' 

_ J2(-l) m+m 'a lm a* m ,Y; m (mrrM)- 



21 + 

mm' 

So the cross-term ai m a* lrrit is responsible for the angular dependency of the parity parameter 
G(l; q). We can also calculate the difference between G(l; q) and g(l) by 

G(l;g)-g(l) X~(l; q) - X+(Z; qj/gg) 

9(1) ~ P-(l) ' 1 } 

From the relation A(/) = 0(±), we know that X ± {l;q) < P~{1), and G( ' ; ^ g(0 < 1 for 
/ > 3. So, we conclude that G{1; q) mainly stands for the amplitude of the original parity 
parameter g(l). At the same time, due to the rotational variance of G(l;q), we can study 
the possible preferred direction, which may reveal hints on the origin of the observed parity 
asymmetry in CMB field. 

Let us show that G(l; q) map depends on the angular q. As we have mentioned, q labels 
the z-axis direction in the rotated coordinate system, and (8, 0) is just the polar coordinate 
of this direction in Galactic coordinate. We plotted the parameter G(l; q) as a function of 
q for 3 < I < 22, and found that G(l; q) have the similar morphology for I > 4, which 
is clearly shown in Fig. [3j (Note that, the morphology of I = 3 map is different, which 
may relate to the unsolved low quadrupole problem as well as the alignment of quadrupole 



1 Throughout this paper, we use the polar coordinate (9, </>) in the Galactic system, which relates to the 
Galactic coordinate (I, b) by I = 90° — 9 and 6 = 0. 
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and octupole ( IBennett et al.ll201ll ).) In Tabled], we list the preferred directions q, where 
the parity parameter G (l\ q) for each I is minimized (note that different from the problem 



in (IGubitosi et al 



20111). here as the wide ly discussed aligment problem of quadrupole and 



octupole in (jde Oliveira-Costa et al.ll2006l ). the uncertainties of the preferred directions are 
difficult to be defined), which are very close with each other for I > 4. So, we can choose the 
special direction q (note that — q is another equivalent preferred direction), where all the 
parameters G(l; q) are minimized or maximized. We picked out these regions, and plotted 
them in the Galactic coordinate system in Fig. |H It is interesting to find that the preferred 
directions q, where parity violation is largest, are coincident with the WMAP7 KD direction 
(IJarosik et al.l 120111 ). while the preferred directions q, where parity asymmetry is smallest, 
are nearly perpendicular to the KD direction. If we assume the parity asymmetry in the 
CMB has the cosmological origin, it is very hard to explain these coincidences. So, the 
coincidence of the preferred direction q with the WMAP7 KD direction implies that the 
CMB parity asymmetry may relate to the possible contamination of residual WMAP KD 
component. 

Although we will not detailedly study the physical mechanism in this paper, we could 
provide some possible explanations for this coincidence problem. It is noticed that there is 
not a great deal of residual dipole in the WMAP data, so a possible explanation could be 
connected to the use of the dipole as a photemetric calibrator in the WMAP data set. 

Another possible explanation is related to the contaminations generated by the collective 
emission of Kuiper Belt Objects (KBOs) and other minor bodies in the solar system where 
the KD direction is localized. Since the emission of KBOs is nearly independent of the 
frequency in the WMAP freq uency range, this contamination is very hard to be removed in 
the WMAP data analysis. In (IMaris et al.l 1201 ll ; lHansen et al.ll2011bl ). it was discussed that 
this foreground residual could leave significant parity asymmetry in the CMB data. 

Besides, the explanation may also relate to the measure deviation of the WMAP kine- 
matic dipole, which could be caused by the measure error in di pole direction, antenna point- 
ing direction, sidelobe pickup contamination, and so on. In (jLiu fc Lill201ll ). it was found 
that this KD deviation could generate the artificial CMB anisotropics in the low multipoles. 
If this is true, these artificial components may account for the CMB parity violation. 

In order to cross- check this result, we consider another rotationally variant estimator, 
which is proposed by (jde Oliveira-Costa et al.ll2006l ) 



D(l) 



21 



1 1 

— y 

+i ^. 



2 i 1 2 

m \aim\ ■ 



(11) 



m=—l 



If our Universe is statistically isotropic, the ensemble average of this estimator is related to 
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the power spectrum as follows: 



W)) = (12) 



As well discussed, this statistic has also chosen the preferred direction, i.e. the z-axis direc- 
tion. In addition, this statistic favors high ms and so it works well in searches for planarity. 
In a quantum mechanical system, this qua ntity also corresponds to the angular momentum 



along the z-axis direction (lEdmondsl 119851 ) . Due to the rotational variance of this quantity, 
we can define D(l; q) and the corresponding parity parameter G(l; q), where q is the Euler 
rotation angle. We notice that in the definition of G(l; q), due to the factor m 2 in D(l), the 
weight of higher multipoles are much greater than the lower ones. So in this paper, we only 
consider the parity parameter G(l; q) for / < 10, where the parity violation is obvious. We 
plot the quantity G(l; q) as a function of q for low multipoles in Fig. and find that these 
maps have the quite similar morphology, especially for 4 < / < 10. The preferred direc- 
tions q (similar, — q is another equivalent preferred direction), where the parity parameter 
G{1] q) for each I is minimized, are also listed in Table [2j Again, we find the similar results: 
for 4 < I < 10, the bluer regions (larger parity violation) are excellently coincident to the 
WMAP7 KD direction, while the redder regions (smaller parity asymmetry) are nearly per- 
pendicular to the KD direction. So this cross-check agrees with our previous finding: the 
CMB parity asymmetry in the low multipoles may connect with the possible contamination 
of WMAP KD component. 



4. Conclusion 

In this paper, we have investigated the directional properties of CMB parity asymmetry. 
In order to break the rotational invariance of the CMB power spectrum, we defined two dif- 
ferent power spectrum estimators so that a special direction is picked out. By rotating these 
estimators with respect to the Galactic coordinate system, we studied the corresponding 
parity parameters as functions of the preferred directions q, where the parity parameters are 
minimized or maximized. We found that these preferred directions are aligned (parallel or 
perpendicular) with the WMAP7 kinematic dipole direction, which implies that the CMB 
parity asymmetry may be produced by the systematics associated with kinematic dipole. 
This study also shows that the effect of the WMAP kinematic dipole may extend to the 
higher multipoles I ~ 22. 

The Planck surveyor possesses wide frequency coverage and systematics distinct from 
the WMAP. In particular, it may take advantage of both COBE and WMAP results for 
the dipole calibration and meanwhile Planck satellite will have high signal-to-noise ratio in 
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polarization data. Therefore, we may apply the similar tests on the CMB TE and EE data 
from the Planck surveyor, and hope to resolve the association between the parity asymmetry 
and kinematic dipole. 
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Fig. 1. — The theoretical (blue curve) and observed (red curve) values of C(9 — tt) as a 
function of l' min , where the WMAP7 power spectra C{1) have been used as the observed 
data. The error bars indicate the la confident levels caused by cosmic variance. 
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Fig. 2. — The TT correlation function estimated from WMAP7 observational data sets 
and a theoretical prediction with la interval (shaded with Cyan color). We estimated the 
theoretical prediction respectively with l' min = 2 and l' min = 1, where we used the residual 
dipole anisotropy of the ILC7 map. The foreground-contaminated region in the data sets is 
excluded by the WMAP KQ75 mask. 



13 



1=3 



1 = 4 




1=7 



1=8 



1=11 



1=12 




1=21 



1 = 22 




Fig. 3. — The parameter G(l; q) based on the estimators in Eq. (J7J) as a function of q 
for 1 = 3, 7, 11, 21 (left panels) and 1 = 4, 8, 12, 22 (right panels). 
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Fig. 4. — The ILC7 dipole component in the Galactic coordinate system. The functions 
G(l; q) (4 < I < 22) minimize at the white regions, and maximize at the black regions. Note 
that the center direction of the white regions are q = (46.59°, 277.98°) and — q, while those 
of the black regions are q = (50.55°, 167.34°) and — q. 
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Fig. 5. — The parameter G(l; q) based on the estimators in Eq. ljTTj) as a function of q 
for I = 3, 5, 7, 9 (left panels) and I = 4, 6, 8, 10 (right panels). 
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Table 1: The WMAP7 kinematic dipole direction is compared with the preferred direction 
q = (0,<p), where the parity parameter G(l; q) (based on the estimator in Eq. ((?])) is 
minimized. Note that — q is another preferred direction. 







4> 



cos a 



KD 

I = 3 
/ = 4 
/ = 7 
/ = 8 
/ = 11 
/ = 12 
I = 21 
I = 22 



41.74 
85.22 
46.59 
48.19 
48.99 
49.77 
49.77 
51.32 
50.50 



263.99 
204.61 
280.89 
279.14 
277.03 
277.73 
277.73 
283.36 
284.06 



0.400 
0.975 
0.976 
0.979 
0.976 
0.976 
0.957 
0.957 



a a is the angle between q and the KD direction. 
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Table 2: The WMAP7 kinematic dipole direction is compared with the preferred direction 
q = (0,<f>), where the parity parameter G(l; q) (based on the estimator in Eq. (II ip ) is 
minimized. Note that — q is another preferred direction. 







4> 



cos a 



KD 


41.74 


263.99 




I = 3 


86.42 


206.02 


0.401 


/ = 4 


45.80 


303.20 


0.890 


I = 5 


48.19 


305.86 


0.867 


/ = 6 


52.08 


274.22 


0.975 


/ = 7 


57.91 


279.84 


0.939 


I = 8 


39.20 


255.57 


0.994 


I = 9 


37.20 


252.90 


0.989 


I = 10 


40.30 


249.17 


0.985 



a a is the angle between q and the KD direction. 



